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Abstract. Several spatially continuous pedestrian dynamics models have 
been validated against empirical data. We try to reproduce the experi- 
mental fundamental diagram (velocity versus density) with simulations. 
In addition to this quantitative criterion, we tried to reproduce stop-and- 
go waves as a qualitative criterion. Stop-and-go waves are a characteristic 
phenomenon for the single file movement. Only one of three investigated 
models satisfies both criteria. 



1 Introduction 

The applications of pedestrians' dynamics range from the safety of large events 
to the planning of towns with a view to pedestrian comfort. Because of the 
computational effort involved with an experimental analysis of the complex col- 
lective system of pedestrians' behavior, computer simulations are run. Models 
continuous in space are one possibility to describe this complex collective system. 

In developing a model, we prefer to start with the simplest case: single lane 
movement. If the model is able to reproduce reality quantitatively and quali- 
tatively for that simple case, it is a good candidate for adaption to complex 
two-dimensional conditions. 

Also in single file movement pedestrians interact in many ways and not all 
factors, which have an effect on their behavior, are known. Therefore, we follow 
three different modeling approaches in this work. All of them underlie diverse 
concepts in the simulation of human behavior. 

This study is a continuation and enlargement of the validation introduced 
in [7]. For validation, we introduce two criteria: On the one hand, the relation 
between velocity and density has to be described correctly. This requirement is 
fulfilled, if the modeled data reproduce the fundamental diagram. On the other 
hand, we are aiming to reproduce the appearance of collective effects. A char- 
acteristic effect for the single file movement are stop-and-go waves as they are 
observed in experiments [S]. We obtained all empirical data from several exper- 
iments of the single file movement. There a corridor with circular guiding was 
built, so that it possessed periodic boundary conditions. The density was var- 
ied by increasing the number of the involved pedestrians. For more information 
about the experimental set-up, see [5], [6]. 



2 Spatially Continuous Models 



The first two models investigated are force based and the dynamics are given by 
the following system of coupled differential equations 

m i d ^=F i with F l = F^ + F[ e P and ^ = v t (1) 
dt dt 

where Fi is the force acting on pedestrian i. The mass is denoted by mi, the 
velocity by V{ and the current position by £j. This approach derives from sociol- 
ogy [1] . Here psychological forces define the movement of humans through their 
living space. This approach is transferred to pedestrians dynamics and so Fi is 
split to a repulsive force F[ ep and a driven force Ff rv . In our case the driving 
force is defined as 

where vf is the desired velocity of a pedestrian and r their reaction time. 

The other model is event driven. A pedestrian can be in different states. A 
change between these states is called event. The calculation of the velocity of 
each pedestrian is straightforward and depends on these states. 



2.1 Social Force Model 

The first spatially continuous model was developed by Helbing and Molnar [3] 
and has often been modified. According to [5] the repulsive force for pedestrian 
i is defined by 

prep^ =J2fi j (x i ,x j ) + Zi(t) with f ij (x i ,x j ) = -d x A(Ax i , j -D)- B , (3) 

with A = 0.2, B = 2,D = 1 [m],r = 0.2 [s] and A the distance between 

pedestrians i and j. The fluctuation parameter £,•(<) represents the noise of the 
system. In two-dimensional scenes, this parameter is used to create jammed 
states and lane formations [2]. In this study, we are predominantly interested 
in the modeled relation between velocity and density for single file movement. 
Therefore the fluctuation parameter has no influence and is ignored. 

First tests of this model indicated that the forces are too strong, leading to 
unrealistically high velocities. Due to this it is necessary to limit the velocity to 
v-max, as it is done in [3] 

l+ \ _ / Vi(t), if \vi(t)\ < 
\vmax, otherwise 

In our simulation we set v max = 1.34 [^]. 



2.2 Model with Foresight 



In this model pedestrians possess a degree of foresight, in addition to the current 
state of a pedestrian at time step t. This approach considers an extrapolation of 
the velocity to time step t + s. For it [8J employs the linear relation between the 
velocity and the distance of a pedestrian i to the one in front Ax^i+ i(t). 

Vi(t) = a Axi t i+i(t) - b (5) 

For a = 0.94 [-y] and b = 0.34 [i] this reproduces the empirical data. So with 
(O Vi(t + s) can be calculated from Axa+i(t + s) which itself is a result of the 
extrapolation of the current state 

Axi, i+ i (t) + Av(t) s = Axi, i+ i(t + s) (6) 

with Av(t) — Vi+\(t) s — Vi(t) s. Finally, the repulsive force is defined as 

F[ ep (t) = _ v i- y i( t + s ) . (7) 

Obviously the impact of the desired velocity v® in the driven force is negated by 
the one in the repulsive term. After some simulation time, the system reaches 
an equilibrium in which all pedestrians walk with the same velocity. In order to 
spread the values and keep the right relation between velocity and density, we 
added a fluctuation parameter uniformly distributed in the interval 

[—20, 20] reproduced the scatter observed in the empirical data. 



2.3 Adaptive Velocity Model 

In this model pedestrians are treated as hard bodies with a diameter di pQ . The 
diameter depends linearly on the current velocity and is equal to the step length 
st in addition to the safety distance f3 

di (t) = e + fvt (t) = sti (t) + ^ (t) . (8) 

Based on [9] the step length is a linear function of the current velocity with 
following parameters: 

sti(t) = 0.235 [m] +0.302 [s]vi{t) . (9) 

e and / can be specified through empirical data and the inverse relation of ([S]). 
Here e is the required space for a stationary pedestrian and / affects the velocity 
term. For e = 0.36 [m] and / = 1.06 [s] the last equations ([8} and (|9|) can be 
summarized to 

Pi(t) = di(t) - sU(t) = 0.125 [to] + 0.758 [s] v z (t) . (10) 
By solving the differential equation 

^ = F drv = — V{t) => v(t)=v + ccxp(--) 1 iovceTR, (11) 
at t \ t J 



(a) Demonstration of the (b) Conception of the adaptive velocity 

parameters d, st and ft 

Fig. 1. Left: connection between the required space d, the step length st and 
the safety distance f3. Right: The adaptive velocity with acceleration until i^eci, 
than deceleration until t acc i, again acceleration until td eC 2 and so on. 

the velocity function is obtained. This is shown in Fig. Q] together with the 
parameters of the pedestrians' movement. 

A pedestrian is accelerating to their desired velocity v® until the distance to 
the pedestrian in front is smaller than the safety distance. From this time on, 
he/she is decelerating until the distance is larger than the safety distance and so 
on. Via Axi } i+i,di and /3, those events could be defined: deceleration (fl"2")) and 
acceleration (fTB")) . 

To ensure good performance for high densities, no events are explicitly cal- 
culated. But in each time step, it is checked whether an event has taken place 
and td ec ,t acc or t co u are set to t accordingly. The time step, At, of 0.05 seconds 
is chosen, so that a reaction time is automatically included. The discrete time 
step could lead to configurations where overlapping occurs. To guarantee volume 
exclusion, case (fT4)) is included, in which the pedestrians are too close to each 
other and have to stop. 
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3 Validation with Empirical Data 

For the comparison of the modeled and experimental data, it is important to 
use the same method of measurement. [S] shows that the results from different 
measurement methods vary considerably. The velocity Uj is calculated by the 
entrance and exit times tf 1 and t° ut to two meter section. 



To avoid discrete values of the density leading to a large scatter, we define the 
density by 

P{t) = ^=f^\ (16) 
2 m 

where &i(t) gives the fraction to which the space between pedestrian i and i + 1 
is inside the measured section, see [6]. pi is the mean value of all p(t), where 
t is in the interval [t\ n ,t° ut ]. We use the same method of measurement for the 
modeled and empirical data. The fundamental diagrams are displayed in Fig. [21 
where N is the number of the pedestrians. 
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(a) Social force model 
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(b) Model with foresight 
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(c) Adaptive velocity model 



(d) Empirical data 



Fig. 2. Validation of the modeled fundamental diagram with the empirical data 
(down right) for the single file movement. 



The velocities of the social force model are independent of the systems density 
and nearly equal to the desired velocity Vi = v® ~ 1.24 [—]. Additionally we 
observe a backward movement of the pedestrians and pair formation. Because 
of these unrealistic phenomena are not observed in the other models, we suggest 
that this is caused by the combination of long-range forces and periodic boundary 
conditions. 

In contrast, the model with foresight results in a fundamental diagram in 
good agreement with the empirical one. Through the fluctuation parameter, the 
values of the velocities and densities vary as in the experimental data. 

We are satisfied with the results of the adaptive velocity model. For reducing 
computing time, we also tested a linear adaptive velocity function, leading to a 
70% decrease in computing time for 10000 pedestrians. The fundamental dia- 
gram for this linear adaptive velocity function is not shown, but also reproduces 
the empirical one. 



4 Reproduction of Stop-and-Go Waves 

During the experiments of the single file movement, we observed stop-and-go 
waves at densities higher than two pedestrians per meter, see Fig. 5 in [5]. There- 
fore, we compare the experimental trajectories with the modeled ones for global 
densities of one, two and three persons per meter. The results are shown in FigEl 
Since the social force model is not able to satisfy the criterion for the right rela- 
tion between velocity and density, we omit this model in this section. Figure [3] 
shows the trajectories for global average densities of one, two and three persons 
per meter. From left to right the data of the model with foresight, the adaptive 
velocity model and the experiment are shown. 

In the experimental data, it is clearly visible that the trajectories get unstead- 
ier with increasing density. At a density of one person per meter pedestrians stop 
for the first time. So a jam is generated. At a density of two persons per me- 
ter stop-and-go waves pass through the whole measurement range. At densities 
greater than three persons per meter pedestrians can hardly move forward. 

For the extraction of the empirical trajectories, the pedestrians' heads were 
marked and tracked. Sometimes, there is a backward movement in the empirical 
trajectories caused by self-dynamic of the pedestrians' heads. This dynamic is 
not modeled and so the other trajectories have no backward movement. This 
has to be accounted for in the comparison. 

By adding the fluctuation parameter Q(t) £ [—20,20] to the model with 
foresight a good agreement with experiment is obtained for densities of one and 
two persons per meter. The irregularities caused by this parameter are equal to 
the irregularities of the pedestrians dynamic. Nevertheless, this does not suffice 
for stopping so that stop-and-go waves appear, Fig |3(d)| and Fig 3(g) 



With the adaptive velocity model stop-and-go waves already arise at a density 
of one pedestrian per meter, something that is not seen in experimental data. 
However, this model characterizes higher densities well. So in comparison with 
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(a) Model with foresight (b) Adaptive velocity model (c) Experimental data 
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(d) Model with foresight (e) Adaptive velocity model (f) Experimental data 
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(g) Model with foresight (h) Adaptive velocity model (i) Experimental data 



Fig. 3. Comparison of modeled and empirical trajectories for the single lane 
movement. The global density of the system is one, two or three persons per 
meter (from top to bottom). 



Fig. 3(h) and Fig. 3(i) the stopping-phase of the modeled data seems to last for 
the same time as in the empirical data. But there are clearly differences in the 
acceleration phase, the adaptive velocity models acceleration is much lower than 
seen in experiment. 

Finally other studies of stop-and-go waves have to be carried out. The occur- 
rence of this phenomena has to be clearly understood for further model modi- 
fications. Therefore it is necessarry to measure e. g. the size of the stop-and-go 
wave at a fixed position. Unfortunately it is not possible to measure over a time 



interval, because the empirical trajectories are only available in a specific range 
of 4 meters. 

5 Conclusion 

The well-known and often used social force model is unable to reproduce the 
fundamental diagram. The model with foresight provides a good quantitative 
reproduction of the fundamental diagram. However, it has to be modified fur- 
ther, so that stop-and-go waves could be generated as well. The model with 
adaptive velocities follows a simple and effective event driven approach. With 
the included reaction time, it is possible to create stop-and-go waves without 
unrealistic phenomena, like overlapping or interpenetrating pedestrians. 

All models are implemented in C and run on a simple PC. They were also 
tested for their computing time in case of large system with upto f 0000 pedes- 
trians. The social force model offers a complexity level of 0(N 2 ), whereas the 
other models only have a level of O(N). For this reason the social force model 
is not qualified for modeling such large systems. Both other models are able to 
do this, where the maximal computing time is one sixth of the simulated time. 

In the future, we plan to include steering of pedestrians. For these models 
more criteria, like the reproduction of flow characteristics at bottlenecks, are 
necessarry. Further we are trying to get a deeper insight into to occurrence of 
stop-and-go waves. 
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